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We study the distributions functions for global partial density of states (GPDOS) in quasi-one- 
dimensional (Q1D) disordered wires as a function of disorder parameter from metal to insulator. We 
consider two different models for disordered Q1D wire: a set of two dimensional S potentials with an 
arbitrary signs and strengths placed randomly, and a tight-binding Hamiltonian with several modes 
and on-site disorder. The Green functions (GF) for two models were calculated analytically and it 
was shown that the poles of GF can be presented as determinant of the rank N x N, where N is 
the number of scatters. We show that the variances of partial GPDOS in the metal to insulator 
crossover regime are crossing. The critical value of disorder w c where we have crossover can be used 
for calculation a localization length in Q1D systems. 

PACS numbers: 72.10.Bg, 72.15.Rn, 05.45.-a 
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I. INTRODUCTION 



Calculation of density of states (DOS) allowed us ob- 
tain many properties of the system under consideration, 
such as charging effects, electrical conduction phenom- 
ena, tunneling spectroscopy or thermodynamic proper- 
ties. Furthermore, the decomposition of DOS in partial 
density of states (PDOS) and global PDOS (GPDOS), 
which appear naturally in scattering problems in which 
one is concerned with the response of the system to small 
perturbation SU(x) of the potential U(x), plays an im- 
portant role in dynamic and nonlinear transport in meso- 
scopic conductorsi2iMi£i£iL&. Particularly the emissiv- 
ity, which is the PDOS in configuration space for elec- 
trons emitted through arbitrary le&c&LHL, always present 
in physical phenomena where quantum interference is im- 
portant. As shown inii the heat flow, the noise prop- 
erties of an adiabatic quantum pump can expressed in 
terms of generalized parametric emissivity matrix v[X] 
(the diagonal element v aa [X] of which is the number of 
electrons entering or leaving the device in response to 
small change SU (x) , such as a distortion of the confining 
potential). The nondiagonal element v a p[X] of paramet- 
ric emissivity matrix determines the correlation between 
current in the contacts a and (3 due to a variation of 
parameter X—. Note that the elements of GPDOS are 
closely related to characteristic times of the scattering 
process, consequently, to the absolute square of the scat- 
tering states. Particularly in ID systems v a p and v aa are 
related to Larmor transmitted time tt (or Wigner delay 
time) and reflected time weighted by the transmission 
coefficient T^*i£ and reflection coefficient R, respectively. 
As it was mentioned by Biittiker and Christen^ the dy- 
namic response of the system to an external time de- 



pendent perturbation, i.e., the emittance in general not 
capacitance-like, i.e. the diagonal and the off-diagonal 
omittance elements are not positive and negative val- 
ues, respectively. Whenever the transmission of carriers 
between two contacts predominates the reflection, the 
associated emittance element changes sing and behaves 
inductance-like. This type of cross over behavior for di- 
agonal element of emittance v aa (taking into account the 
Coulomb interaction of electrons inside the sample) was 
found inia where they study the distribution function 
(DF) of emittance. They have found that in the range 
of weak disorder, when the system is still conducive the 
DF is Gaussian-like. With increasing disorder the DF 
becomes non-Gaussian. 

The purpose of this paper is to study numerically the 
behaviors of DF of diagonal v aa and off-diagonal ele- 
ments v a p of global PDOS in the Q1D disordered wires, 
where not so much known about the DF. We study 
three different regimes of transport: metallic (£ >> L), 
where £ is the localization length and L the typical 
size of the system, insulating (£ << L) and crossover 
(£ *~ L). We show that in intermediate regime of trans- 
port between the metallic and insulating regimes there 
is the critical value of disorder w c when we observe 
cross over between the variances var(^ Qa ) and \&i:{v a p) 
(see Fig.l). This critical w c determines the localization 
length of Q1D system for given length L and number of 
modes M. It turns out that in metallic regime P{v a p) 
is Gaussian which means that the first and the second 
moments (i.e., the average < v a p > and the variance 
var(i/ Q( g) =< > — < v a p > 2 ) are enough to describe 
the behavior of P(v a p)- I n the strong localization regime 
the distribution of v a p is log normal, which means that 
the lni/ Q/ 3 follows a Gaussian distribution. As regards 
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the distribution function of v aa we can say that in the 
strong localization regime it characterized by an expo- 
nential tail, the values of positive and that the 
dynamic response of the system is capacitivelikei*. In 
the metallic regime the emittance has non Gaussian-like 
behavior and some of the value of v aa are negative (in- 
ductivelike behavior 

The paper is organized as follows. In the next section 
we present our model and set the basis on numerical cal- 
culation for obtaining the probability distributions of v aa 
and v a j3 for different regimes. In section Hill we study the 
behavior of var(V QQ ,) and vai(y a p) as a function of disor- 
der strength w. In section llVl we calculate the distribu- 
tion functions for v aa and v a p in three different regimes 
of transport mentioned in Introduction. The paper is 
included in section Ivl 



with V n , x n and y n be arbitrary parameters and Q1D 
lattice of size L x W (L 3> W, L is the length and W is the 
width of the system), where the site energy can be chosen 
randomly. In both cases analytically we have calculated 
the Green's function of 1QD (— ) and use them in our 
numerical calculations (see Appendix). The elements of 
global PDOS u a p, in the case of a tight-binding model 
can be calculate in terms of the scattering matrix and 
the Green Function. To calculate the scattering matrix 
elements, corresponding to transmission between modes 
n and m, we start from the Fisher-Lee relationi^i^, which 
expresses these elements in terms of the Green's function: 



II. THE MODEL AND NUMERICAL 
PROCEDURES 

The localization length £ is obtained from the decay of 
the average of the logarithm of the conductance, In g, as 
a function of the system size L 



£ 1 = - lim -J- < hi,? > 

L — >oo LLi 



(i) 



+ iHy/v n v m ^ Xn (ro, ) G (r 0l , r 0j ) \ m (r 0] ) . 

(4) 

Xm \T0j ) is the transverse wavefunction corresponding to 
mode m at the point r Qj and G (r 0i ,r 0j ) is the Green's 
function (GF) for non coinciding coordinates. v m is 
the velocity associated with propagating mode m. The 
LPODS is directly connected to the S-matrix elements 
s nm through the expression^: 



where g is given by Biittiker-Landauer formula—^ (T nm 
is the transmission coefficient from mode n to mode m) 



f3 j / j •*■ nm i 

n.m 



h 



(2) 



We will consider two models: Q1D wire with the set of 5 
scattering potentials of the form 



JV 



dn nm (r) 
dE 



Atti 



Ss r , 



'SU(r) SU(r) 



(5) 



V(x, y) = ^ V n S(x — x n )S(y — y n ), (3) Insertion of Eq. Q in Eq. ©gives: 



dn r . 



dE 



-(r) 



4tt 



i,3 



iXn (rOi)G(r 0i ,r)G(r, r- 0j )x m (r 0j ) + H.c.) 



(6) 



r 



where H.c. denotes Hermitian conjugate. To arrive 
the above expression we have calculated the functional 
derivative of the Green's function by adding to the Hamil- 
tonian of our system the local potential variation SU (r) — 
5U a S(r — r a ) (5U a — > 0), which lead us to the relation^ 



cles of our system: 

dN n 



E 



dn nm (r k ) 



dE V dE 

k 



(7) 



SG(r n ,rr, 
6U(r) 



G(r n ,r)G(r, r rn ). 



After summation over the indices i,j and r k the above 
equation in matrix form can be presented: 



Once we have calculated the local PDOS we can obtain 
the global PDOS adding the local PDOS over the parti- 
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where Q n m matrix defined as 



M 



— Xn J ^ GmjG 
J =1 




(9) 



Xn is the column matrix: 



/ Xn(l) 



V Xn(M) 



(10) 



Here Xm 1S the transpose of the column matrix Xm and 
G m j is the matrix of M x M rank (M is the number of 
modes in each lead, see Appendix). 

Finally, v aa and v a /3 one can get from global PDOS 
dN nm /dE by summing every mode n in lead a and every 
mode m in lead [3 respectively: 



E 



<i/V„ 



dE 



'a/3 



E 



diV„ 



dE 



(11) 



(12) 



Similarly can be written also vpp and vp a and so global 
DOS v must be sum of all GPDOS: 



(13) 



In the case of the Q1D wire with the set of S potentials 
(see Eq. m quantities Q), ifTTj) and calcu- 

lated for tight-binding model one must replace the sign 
of summation by appropriate spatial integration. 

For numerical study we consider a quasi-one dimen- 
sional lattice of size L x If (L » W), where L is the 
length and W is the width of the system. The standard 
tight-binding Hamiltonian with nearest-neighbor interac- 
tion 



H 



e t | n >< n \-t^2\ n >< r 

i,3 



(14) 



where ej is the energy of the site i chosen randomly be- 
tween (— 3f,Tf) with uniform probability. The double 
sum runs over nearest neighbors. The hopping matrix 
element t is taken equal to 1, which sets the energy scale, 
and the lattice constant equal to 1, setting the length 
scale. The energies are measured with respect to the 
center of the band so we will always deal with propa- 
gating modes. Finally our sample is connected to two 
semi-infinite, multi-modes leads to the left (lead a) and 
to the right (lead (3). For simplicity we take the numbers 
of modes in the left and right leads to be the same (M) 
and so the width of this system w becomes equal M (for 
a tight-binding model the numbers of modes coincides 
with the number of sites in the transverse direction). The 



conductance of a finite size sample depends on the prop- 
erties of the system and also on the leads which must be 
taken into account in a appropriate way. In order to take 
into account the interaction of the conductor with the 
leads we introduce a self-energy term " A" as an effective 
Hamiltonian, which will be calculated as (see, e.g>i£) 



ik n a 



Xn(r 0k ) 



A — Ap + A q 



(15) 



(16) 



(17) 



Finally for numerical calculation of DF of vax(y aa ) and 
var(i/ Q/ 3) for M ^> 1 and for higher dimension of the 
system we calculate the Green function as: 



G= [EI-H-A]- 1 



(18) 



To perform numerical calculation of the elements of this 
Green's matrix we will use Dyson's equation, as ini2*2S, 
propagating strip by strip. This drastically reduced 
the computational time, because instead to inverting an 
I? x M 2 matrix, we just have to invert L times L x M 
matrices. In this way we build the complete lattice start- 
ing from a single strip, and introducing one by one the 
interaction with the next strip. Each time we introduce 
a new strip we apply the recursion relations of Dyson's 
equation, until wc finally obtain the Green function for 
the complete lattice. Once we have the Green's function 
matrix we calculate \ax{y aa ) and var(z/ a ^) according to 
Eqs. I|ll|) an d l|12|l. and obtain their probability distri- 
butions for random potentials. Over 250000 independent 
impurity configurations where averaged for each N . 



III. VAR(^ a ) AND VAR(!/ a(S ) VS w 

In this section we are going to study the dependance 
of the var(f Qa ) and var(^ a ^) vs disorder w and vs the 
number of mode M. In Fig.l we show the behavior of 
var(^ QQ ) and var(i/ Q/ 3) as a function of the disorder w. 
Plot is for a sample of L = 400 and M — 4. The crossover 
define a critical value of the disorder w c . In Fig. 2 we show 
the dependence of the critical value w c with the number 
of modes M for several samples. As one can see with 
increasing the number of modes the crossing point moves 
to the left and the w c decreases. This means that in 
the weak localized regime, in analogy with ID systems 
the ratio of localization length £ to the longitudinal size 
of the sample L for given modes M follows, in a good 
approximation, a law of the form 



C(M,w c ,E) 



(19) 
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FIG. 1: Crossover between var(z/ aa ) and var(i/ a /3) for sample 
of L = 400 and M = 4. 



FIG. 3: Universal curve w c for several samples as a function 
of 1/L e ff. o is for ID sample; •, ▲ T and * are for quasi ID 
samples with the numbers of mode M = 2, 3, 4, 5, respectively. 
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FIG. 2: Critical value w c for several samples as a function 
of 1/L: o is for ID sample; •, A, T and ★ are for quasi ID 
samples with the numbers of mode M = 2, 3, 4, 5, respectively. 

where C is a constant that depends from the M, w c and 
energy. With appropriate choice of an effective length 
L eff = L(a + bM c ) (with a = 0.967, b = 0.035 and 
c = 2.33) we were able to show that all the curves pre- 
sented in the Fig. 2 collapse into universal curve in Q1D 
system, supporting the applicability of the hypothesis 
of single-parameter scaling^! in disordered systems. In 
Fig. 3 we plot this curve for w c as a function of 1/L e //. 
The different values of modes are specified inside the fig- 
ure. 

In strictly ID system, following^ one can write (i>t = 
v a f3 + Vf} a and v R = v aa + vpp) 

(hxv R ) = (lau) + (hiR) (20) 



(IntAr) = (In j/) + (InT) (21) 

where R and T are the reflection and transmission coef- 
ficients respectively and (...) denotes averaging over en- 
semble. Using the asymptotic behavior of < InT > and 



FIG. 4: Average of Inun and In it as a function L/£. Solid 
curves are given by Eqs. 11221 and 12311 . The data points (•) 
are the numerical results for a sample of L — 400. 



< In R > as L — > oo (see, e.gM-) these expressions in 
weak disorder regime can be rewritten as: 

(In v R ) = (In v) + ln(l - e~ 2L /«) (22) 



(lni*r) = (In i/) -2L/£ (23) 

In Fig. 4 we plot average of h\v R and h\VT for different 
values of disorder w c as a function L/£. We see that nu- 
merical data for these quantities very well coincide with 
Eqs. I|22|) and (|23|l for small w c or for large L/£. 

IV. PLOTS AND DISCUSSIONS 

We are analyzed the DF P(v aa ) and P(v a p) along the 
transition from the metallic to the insulating regime for 
several samples sizes. We found that the relative shape of 
the DF depends only from the disorder parameter L/£, 
i.e. with increasing of the number of modes M we al- 
ways can find an appropriate range of w for which all the 
curves have the same form. Therefore in the rest of the 



5 



section, without losing generality we present our results 
for a sample of L = 400 and M — 4 for several values of 
the disorder w. 

In the metallic regime when the system size much 
smaller than the localization length L << £ the dis- 
tribution functions are shown in Fig. 5 with (W = 0.2, 
Lji = 0.17 and (g) = 2.52). We have checked that 
the distribution of P{v a p) is Gaussian-like and can be 
fit with the following expression (B = 1.0, /i = 116.5 and 
a = 20.2): 



with: 



P{y a p) 



B 



-(v af s-nf/2(T 2 



2vrcr 



(24) 



In spite of the fact that in our numerical studies we deal 
with Q1D systems where the numbers of modes M > 1, 
still the Gaussian-like behavior of the v a p in ballistic 
regime can be understood well if we recall the fact that 
v> a p connected with physically meaning full times charac- 
terizing the tunneling process 5 ". In ID systems GPDOS is 
related to Larmor transmitted time tt (or Wigner delay 
time) weighted by the transmission coefficient 5 ' 12 , 



T 



v a f3 



-t t , 



(25) 



The quantity tt, which links to the density of states of 
the system^ 3 , and can be presented 



tt = him / G(x, x) dx — him ■. 

Jo I "E 



jOlnt 



4E 



(26) 

where G(x, x) is the GF for the whole system, t and r are 
the transmission and reflection amplitudes from the finite 
system, r' is the reflection amplitude of the electron from 
the whole system, when it falls in from the right. 

The second term in Eq. (|26|l becomes important for 
low energies and/or short systems. This term can be 
neglected in the semiclassical WKB case and, of course, 
when r (and so r 1 ) is negligible, e.g., in the resonant 
case, when the influence of the boundaries is negligible. 
Of course the distribution function of i/ a/ g (Eq. is 
affected by correlations between the value of the DOS 
(or Wigner delay time) and the transmission coefficient 
of resonances via localized states, but still it can capture 
some general behavior Wigner delay time in ID system 
in the regime where T ~ 1. Wigner delay time in ID 
and in the ballistic regime is given by Gaussian function 
and can be characterized by a first moment and a second 
cumulant 24 i 25 . 

Similar relation to Eq. (|25|l holds for v aa : 



R 

Vaa — ~~^ T Ri 



(27) 



where tr characterize the reflection time and defined as: 



tr 



= him 



d In r 1 — r 2 — t 2 



dE 



AEr 



(28) 



9(x) 



exp 



dx 



o 2G(x,x) 



We note that for an arbitrary symmetric potential, 
V((L/2) + x) = V((L/2) — x), the total phases accumu- 
lated in a transmission and in a reflection event are the 
same and so the characteristic times for transmission and 
reflection corresponding to the direction of propagation 
are equal 



tt = tr 



(29) 



as it immediately follows from Eqs. Il2(i[l and J2HJ- For 
the special case of a rectangular barrier, Eq. (|29|l was 
first found in RefiSi. Comparison of the Eqs. I|26|) and 
B28J1 shows that for an asymmetric barrier Eq. I)29|l breaks 
dowrtSi. 

As one can see from Fig. 5 in the same regime DF 
P{v a a) include big range of negative v aa values indicat- 
ing a predominantly inductive dynamic response of the 
system to an external ac electric field according^. For 
positive values of v aa the tail of the distribution P{v aa ) is 
fairly log-normal with following parameters (B = 0.875, 
fx = 60.5 and a = 0.25): 



P(v aa ) 



B 



-(lni/ aa -/i) 2 /2a 2 



2ll(TV aa 



(30) 



With increasing the disorder w, when we almost are 
in crossover regime we obtain a wide range variety of 
broad distributions as shown in Fig. 6 where we plot 
DF for two values of disorder: w — 0.5 (L/£ = 0.69 and 
(g) = 0.75) in the left panel and w = 0.6 (L/£ = 0.93 
and (g) = 0.5) in the right panel. As one can see from 
Fig. 6 (right panel) P{v a p) has a flat part for almost in 
all the range of iy a p while in the left panel it has a strong 
decay. In both cases the distributions for P{v a p) can 
be fitted to two log-normal tails. This type of behavior 
is typical also for distribution of conductance g in the 
same range of parameters in Q1D, as one can see from 
the same Fig. 6 where we present P(g). For values g < 
1 we have a plat part while in the regime g > 1 we 
get for distribution of conductance a strong decay. This 
is in complete agreement with a number of numerical 
simulation in the intermediate regime (see e.gm^^^). 

As for P{vaa) it is shifted to right, to much larger 
value of v a a! which means that it becomes less conduc- 
tive. For this range of parameters DF is still quite sym- 
metric (right panel) but more wider if we compare with 
the DF from the Fig. 5. The P(v aa ) for w = 0.6 becomes 
less symmetric (left panel). 

Further increase of the disorder w (in the insulating 
region) P(v a p) becomes a one-side log-normal distribu- 
tion. This type of behavior was predicted for distribu- 
tion of conductance g insSisi and numerically calculated 

j n 29.30.32 

With regard to P(v aa ), we can mention that the tail 
of the distribution follows a power-law decay P{v aa ) oc 
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FIG. 5: Probability distributions of v aa and v a g in the metal- 
lic regime ({g) — 2.52) for a disorder of w = 0.2. The solid 
lines correspond to a gaussian distribution for u a/ s and a log- 
normal tail distribution for 



l/f™ a , with to ~ 2.3. On the other hand, as w increases 
P(va/s) shows a tail in the negative region of f a /3- In 
Fig. 7 we plot the distribution P{y a a) and P(y a j3) for a 
disorder w = 1 (Z/£ = 2.6 and (g) = 0.08). 

Deeply in the localized regime (L ^> £ and (g) « 0) the 
distribution of ^ a /3 is log-normal as one can from Fig. 8 
where we fit P(hn> a p) to a Gaussian distribution: 

P{X) = ^L- e -(ln*-p) 2 /2- 2 ; (31) 

with B = 0.997, v = -460.5 and a = 27.7. 

The shape of P(y a a) is highly asymmetric with two 
peaks very closed each other. The position and the high 
of this peaks depend on the disorder parameter and cause 
several shapes of the distribution function. The tail of the 
distribution follows a power-law decay P(v aa ) cx l/v™ a , 
with to ~ 2.0 



V. CONCLUSIONS 

We study the distributions functions for global par- 
tial density of states in quasi-one-dimensional disordered 



wires as a function of disorder parameter from metal to 
insulator. We consider two different models for disor- 
dered Q1D wire: a set of two dimensional S potentials 
with an arbitrary signs and strengths placed randomly, 
and a tight-binding Hamiltonian with several modes and 
on-site disorder. It was shown that the poles of Green 
functions can be presented as determinant of the rank 
N x N, where N is the number of scatters. We show that 
the variances of partial global partial density of states in 
the metal to insulator crossover regime are crossing. The 
critical value of disorder w c where we have crossover can 
be used for calculation a localization length in Q1D sys- 
tems. With increasing the numbers of mode the crossing 
point moves to the left and the w c decreases. 

In the metallic regime when the system size much 
smaller than the localization length L << £ the distri- 
bution function for P(y a p) is Gaussian-like. In the same 
regime the distribution function of P(y aa ) is include big 
range of negative v aa values indicating a predominantly 
inductive dynamic response of the system to an external 
ac electric field according^. For positive values of v aa 
the tail of the distribution P(y aa ) is fairly log-normal. 
Almost in crossover regime the distribution function for 
P{v a p) can be fitted to two log-normal tails. As for 
P{vaa) it is shifted to right, to much larger value of v aa , 
which means that it becomes less conductive. Further in- 
crease of the disorder w (in the insulating region) P(v a p) 
becomes a one-side log-normal distribution. With regard 
to P{v aa ), we can mention that the tail of the distri- 
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bution follows a power-law decay P(v aa ) cx l/v™ a , with 
to ~ 2.3. Deeply in the localized regime (L ^> £ and 
(g) w 0) the distribution of v a p is log-normal and while 
the shape of P{v a a) is highly asymmetric with two peaks 
very closed each other. The position and the high of this 
peaks depend on the disorder parameter and cause sev- 
eral shapes of the distribution function. 
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FIG. 7: Distributions P(v a a) and P(ln^ a/ 3) in the insulating 
regime ((g) = 0.08) for a disorder of W = 1. We have a 
one-side log-normal distribution for v a p and power-law tail 
for v aa , P{v aa ) oc \jv™ a , with m ~ 2.3 
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VII. APPENDIX: DYSON EQUATION IN Q1D DISORDERED SYSTEM AND THE POLES OF 

GREEN'S FUNCTION 

We consider the Q1D wire with the impurities potential of the form: 

JV 



V(x, V)='J2 V n S(x - x n )6(y - y n ), 



(32) 



where V n , x n and y n are arbitrary parameters. The equation for the Green function with above potential V(x,y) is: 

h 2 ( d 2 d 2 



E 



, , , +V c (y) + V(x,y) 
2m \dx dy* 



G(xy; x'y') = S(x - x')5{y - y'), 



(33) 



where the confinement potential V c (y) depends only on the transverse direction y. The Dyson equation for a Q1D 
wire can be written in the forrn^: 



G ac (x,x')=G Q a (x,x')5 ac + J2 / G° a (x,x")S ab V bd (x")G dc (x",x')dx", 

b,d J 

The matrix V a b{x) elements of the defect potential are: 

r N 

V ab (x) = / X* a (y)V(x, y)x b (y) dy = J2 V ab^ - *n), 



(34) 



(35) 
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FIG. 8: Distributions P(v aa ) and P(ln^ a/ 3) in the insulating 
regime ((g) ~ 0) for a disorder of W = 12. The solid lines 
correspond to a gaussian distribution for In v a p that point out 
a log-normal distribution for v a p. The tail of the distribution 
follows a power-law decay P(u aa ) oc l/v^a, with m ~ 2.0 



and denned as: 



Vab(x) = X* a {yn)V n Xb{Vn) 



(36) 



Details on the calculation of the GF G nm (x, x') of Dyson equation ll.'j.'jll for this case, based on the method devel- 
oped in3ii2& will be done elsewhere^. Here we present main results of calculation which will be used in numerical 
calculations. The pole of GF can be rewritten as a determinant of the rank (MN x MN) (M is the number of modes 
and N is the number of delta potentials) The matrix elements of determinant's (Dmn)iu & r e : 



(DmnU = -16* + (I- S nl ){X nl }{r^}. 



(37) 



Here: 



I = 



(38) 







1 



is unit matrix. The Ith scattering matrix {r^} and {A n ;} matrix are matrices M x M and defined in the following 
way: 



(0 
it 



\r {l) ••• 

x 'Ml 



(0 

1M 



MM 



(39) 
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(i) 



{Xnl} 



The quantities r~mm and r^ m are: 



\ ••• A„; 



.(0 



(M) J \ g ... gifcM|a;„-2;i| 



(40) 



i-E^te(*i,*i) 



(41) 



km 



(42) 



respectively, rlii, (to = 1, 2, ...A/) is the complex amplitude of the reflection of an electron from the isolate potential 
Vi with the coordinates xi,yi. Electron incidents from the normal mode to on the left (right) and reflected normal 



mode to on the left (right). 



(i) 

r 

km 



is the complex reflection amplitude of an electron from the same Vi but it incidents 



from the normal mode to on the left (right) side and reflected normal mode k on the same side: By permuting indexes 
k and m in (|42() one can find the complex amplitude r^l . Note that determinant of the matrix is zero, i.e. 



det{r {l) } = 



(43) 



This is follows from the fact that 



r»„r2-r«>.r2.=0, 



mm kk 



mk km 



which can be checked directly if one used the definitions of l|42|) and l|41|) . The rank (MN x M N) of the above 
determinant (see (|37|) ") . after some mathematical manipulation can be reduced to the determinant of the rank (N x N), 
as in the case of ID chain of arbitrary arranged potentials", with the following matrix elements: 



(D N ) n i = -5 n l + (1 - 5. 



nl ) 



w A « + _L yvV p) A (p) 

11 A nl * m / y ' lp 'pi A nl 
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(44) 



Once we know the explicit form of (Dn) h i, we can calculate the scattering matrix elements without determining the 

exact electron wave function in disordered Q1D wire. For example the transmission amplitude from the set of 

N delta potentials is: 



rp(N) ^ik\ \x N — n (Dn) III 

(Dn)uI 



(45) 



where the matrix {D^ni is obtained from the matrix (-Dat)„z (Eq. (|44JH by augmenting it on the left and on the top 
in the following way: 



(DN)nl — 



-iki \xn — xi I 



„( N > p iki\x N -X! | 



(Dn)uI 



(46) 



The reflection amplitude of electrons from the same set of TV delta potentials is given by: 



R (N) _ {DN)nl 
' ' {D]Sl)nl 



(47) 
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where the matrix (-Djv)ni is obtained from the matrix (-Djv) n ; (Eq. (|44[) ) by augmenting it on the left and on the top: 



(Dn)tiI — 



a ik 1 \x N -x 1 \ 



11 



( N ) ikt\x N -xi\ 



' 11 



(Dn) 



111 



(48) 



It can checked directly that Eq. (|45|l for the case of two point scatterers (i.e. N — 2) and for two modes (M = 2) 
lead us to (ai — X2 — x\) 



(1 + r W)(i + rg) + rgrge^-^ - rgrge^ - e^+^rgrff 



1 T ,( 1 )„( 2 ) p 2i/ciai 1 ,( 1 ) T .( 2 )„2ifc2ai / T ,( 1 ) T .( 2 ) i T ,( 1 ) T .( 2 ) N \ p i(fc 
1 '11 '11 e '22 '22 e I '12 '21 + '21 '12 J e 



i+fc 2 )ai 



which, after appropriate notation used in^i, will coincides with their expression of Tn calculated by transfer matrix 
method. 

For R^P from Eq. |g7J we will get 



JX) _L r ( 2 )/"1 , o r ( 1 )^„2 l fc 1 a 1 , f U) (2) I „( 1 )„( 2 )A p i(fc 1 +fc 2 )ai 
r ll + r ll l 1 T Zr ll J e + ^ r 12 r 21 T r 21 r 12 J e 

2 a t _ (r^r^ +r i 2 \ ) r { ^ e ^+ k ^ 



( 1 ) r .( 2 )„2ifciai _ A 1 ) r ( 2 ) „2ik 
' 22 ' 22 e 



1 — r u 'r\i e 



To close this section let us note that to get the expressions for the pole of the GF Eq. (|37|) . for transmission 
amplitude T^i and for Rii in tight-binding model one must to replace the unperturbed GF for normal mode m 



G m {x, x 1 ) 



-i . 2 , exp \ika\x - x | 



with 



k — 

by the appropriate GF for tight-binding model2&: 

G° m (l,n) = - 

Here (x = (E-e)/B) 



2m (E- E„ 



h 2 



^B 2 -(E-e) 



,\l-n\0 



= hi(x — %\J 1 — x 2 ) 



(49) 



(50) 



(51) 



(52) 



and symbol vl — x 2 denotes the positive square roots. 
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